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Nonlinear behavior of electro-osmosis in dilute non-adsorbing polymer solutions with low salinity is investi¬ 
gated with Brownian dynamics simulations and a kinetic theory. In the Brownian simulations, the hydrody¬ 
namic interaction between the polymers and a no-slip wall is considered with Rotne-Prager approximation of 
Blake tensor. In a plug flow under a sufficiently strong applied electric field, the polymer migrates toward the 
bulk, forming a depletion layer thicker than the equilibrium one. Consequently, the electro-osmotic mobility 
increases nonlinearly with the electric field and gets saturated. This nonlinear mobility qualitatively does not 
depend on the details of rheological properties of the polymer solution. Analytical calculation of the kinetic 
theory for the same system reproduces quantitatively well the results of the Brownian dynamics simulation. 


I. INTRODUCTION 

Electro-osmosis is observed widely in many systems 
such as colloids, porous materials and biomenrbranes. It 
characterizes the properties of interfaces between solids 
and electrolyte solutionsRecently, there have been 
growing interests in applications of electro-osmosis. For 
instance, it is used to pump fluids in microfluidic devices 
since it is more efficient than pressure-driven flow^ Ap¬ 
plication to an electrical power conversion is also very 
fascinating in chemical engineering^ When the elec- 
trokinetic properties of a surface are characterized by 
zeta potential, the Smoluchowski equation is often em¬ 
ployed with measurements of the electro-osmotic or elec¬ 
trophoretic mobilities. However, one has to consider 
the validity of the mentioned equation more seriously. 
It is derived from the Poisson-Boltznrann equation and 
the Newton’s constitutive equation for viscous fluids. 
The zeta potential is defined as the electrostatic po¬ 
tential at the plane where a no-slip boundary condi¬ 
tion is assumed. When these equations are not validat- 
able, the Smoluchowski equation is also questionable. In 
the cases of strong-coupling double layer^ inhomogene¬ 
ity of viscosity and dielectric constant near interfaced 
and non-Newtonian fluids^r—, for example, the Poisson- 
Boltzmann and/or the simple hydrodynamic equations 
sometimes do not work well. 

In order to control the electrokinetic properties of 
charged capillaries, the structures of liquid interfaces con¬ 
tacted with charged surfaces are modified by grafting or 
adding polymers*^ In capillary electrophoresis, for exam¬ 
ple, the electro-osmotic flow is reduced by grafted poly¬ 
mers on the interfaces. About surfaces with end-grafted 
charged and uncharged polymers, several studies have 
been also reportedJ^AS Under a weak applied electric 
field, the grafted polymer still remains in the equilibrium 
configuration and the resultant electro-osmotic velocity 
behaves linearly with respect to the electric fields. To 
measure the mobility of such a surface, hydrodynamic 
screening and anomalous charge distributions due to the 
grafted polymers are important. 14-19 When a sufficiently 


large electric field is applied, the polymers are deformed 
by the flow and electric field, and thus, the electro- 
osmotic velocity becomes nonlinear It should be noted 
that the end-grafted polymers cannot migrate toward the 
bulk since one of the ends is fixed on the surfaces. 

When we add polymers into solutions, a depletion or 
adsorption layer is often formed near a solid wall as well 
as diffusive layers of ions in equilibrium states. The in¬ 
teraction between the polymers and the wall determines 
whether the polymers are depleted from or adhere to the 
surfaces. The thickness of the depletion or adsorption 
layers is of the same order of the gyration length of poly¬ 
mers. When polymers adhere to the wall, the viscosity 
near the wall becomes large, so that the electro-osmosis 
mobility is much suppressed^ Moreover, it is known that 
an adsorption layer of charged polymers can change the 
sign of the mobilityiliSlr— The curvature of the surface 
also modulates the surface charge density, and even in¬ 
creases the mobility beyond the suppression caused by 
the viscosity enhancement^ 

Electro-osmosis of a non-adsorbing polymer solution 
was analyzed by two length scales; the equilibrium de¬ 
pletion length S 0 and the Debye length In the de¬ 

pletion layer, the viscosity is estimated approximately by 
that of the pure solvent and it is smaller than the solution 
viscosity in the bulk. When the Debye length is smaller 
than the depletion length, the electro-osmotic mobility is 
larger than that estimated by the bulk value of viscosity. 
Typically for 10 mM electrolyte solutions, one has A ss 
3nm and <5o ~ 100 nm. In such a case^Sr— an electro- 
osmotic flow of high shear rate is localized in the distance 
A from the wall. Thus, the electro-osmotic flow profile 
and resultant electro-osmotic mobility are almost inde¬ 
pendent of the polymers. Actually, such behaviors are 
experimentally observed in solutions of carboxymethyl 
cellulose with urea^27 On the other hand, in the solu¬ 
tions of small polymers with low salinity, typically for 
0.1 mM electrolyte solutions A ~ 30 nm and <5o ~ 5 nm, 
the electro-osmotic mobility is suppressed by the poly¬ 
meric stress 

When a sufficiently strong electric field is applied, the 
electro-osmosis of a polymer solution shows non-linear 
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behaviors j 27,29 These nonlinearities are theoretically an¬ 
alyzed by models of uniform non-Newtonian shear thin¬ 
ning fluids. 9-12 Assuming that polymers still remain lo¬ 
calized in interfacial layers and the viscosity depends on 
the local shear rate as in power-law fluids, their phe¬ 
nomenological parameters are different from those in the 
bulk since the concentration in the interfacial layers is 
different from the bulk concentration. 27 Thus, under¬ 
standing of nonlinear electro-osmosis still remains phe¬ 
nomenological. Furthermore, when shear flow is applied 
to polymer solutions near a wall, it is experimentally and 
theoretically confirmed that cross-stream migration is in¬ 
duced toward the bulker— The concentration profiles 
of the polymer near the wall are calculated and the de¬ 
pletion length dynamically grows tenfold larger than the 
gyration radius^i However, these hydrodynamic effects 
in electrokinetics have not been studied so far to the best 
of our knowledge. 

In this context, the present paper discusses another 
origin of nonlinearity which is induced by the hydrody¬ 
namic interaction between the polymer and wall. Mainly 
the situation of do A is concerned. For this purpose, 
this paper is organized as follows. Section [TT] presents 
a toy model for nonlinear electro-osmosis of dilute poly¬ 
mer solutions. Section EH describes Brownian dynamics 
simulation. Section IIVI presents results of the simula¬ 
tion. Section 0 discusses analytical approach for non¬ 
linear electro-osmosis by using kinetic theory of cross¬ 
stream migration. 31 - Section fVTl outlines the main conclu¬ 
sions. 


II. A TOY MODEL 


First, we propose a toy model for electro-osmosis of 
polymer solutions. A dilute solution of non-adsorbing 
short polymers is considered. The viscosity of the solu¬ 
tion is given by 


V = ?7o(l +%>), 


(1) 


where r/o is the viscosity of the pure solvent, and rj sp is the 
specific viscosity of the solution. The gyration length of 
the polymers is defined as <5o, which is of the same order of 
the equilibrium depletion length. It is assumed that the 
short polymers have ~ 100 nm. Ions are also dissolved 
in the solution with the Debye length A. When a well 
deionized water is considered, the Debye length is of the 
order of A « 10 3 nm although such a salt-free water is 
hardly realized owing spontaneous dissolutions of carbon 
dioxides. The interfacial structure near a charged surface 
is characterized by A and Sq. When an external electric 
field is applied, a shear flow is locally imposed within the 
distance A from the wall, and the resultant shear rate is 


7 ~ 


Mo E 
A 


( 2 ) 


where /io is the electro-osmotic mobility for the pure sol¬ 
vent and is estimated typically as /x 0 ~ 10 -8 m 2 /(V- s). 




FIG. 1. (a) The depletion length as a function of the applied 

electric field, (b) The electro-osmotic mobility as a function 
of the applied electric field. 


According to the studies of the cross-stream migration 
in the uniform shear flow^i, the depletion layer thickness 
depends on the shear rate, 

S « JoOi) 2 , (3) 

where r is the characteristic relaxation time of the poly¬ 
mers, 


VqSq 

k B T ’ 


and is typically 10 4 s. Using eqs. m and ©, the deple¬ 
tion length in the presence of the applied electric field E 
can be expressed by 



for E < Eg, 
for E 0 < E < Ei, 
for E\ < E, 


( 5 ) 


where Eq = X/t^lq, and Ei = Eq^/X/Sq- Here, for sim¬ 
plicity, we assume that the depletion length does not ex¬ 
ceed the Debye length. The effective viscosity in the dou¬ 
ble layer is given by 


VeS ~ Vo 


1 + Vsp ( 1 ^ 


and the nonlinear mobility can be estimated by 
Vo(vo/Ves)- Therefore, the mobility is obtained as 


Vo 


1 + 7?sp(l — (V A)) 

_Mo_ 

l + 7? sp (l ~{E/EiY) 


for E < E 0 , 
for Eq < E < Ei, 


( 6 ) 
M « 


( 7 ) 


Mo for Ei < E. 

Fig- HI (a) shows schematically the thickness of a deple¬ 
tion layer as a function of electric field strength. Fig. [1] 























a dumbbell 


satisfies the fluctuation-dissipation relation as 

(tnaMUfiP)) = 2k B TG%p(t - t'). 
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(a) (b) 




(9) 

To include the effects of the no-slip boundary, Rotne- 
Prager approximation for Blake tensor— is used for the 
mobility tensor for distinct particles (n ^ ro) , 35 i 36 al¬ 
though it is valid only for particles separated far away. In 
this study, we neglect lubrication corrections for nearby 
particles. 37 The Blake tensor for the velocity at X 2 in¬ 
duced by a point force at X\ with the no-slip boundary 
at 2 = 0 is given by the Oseen tensor and the coupling 
fluid-wall tensor as^ 

G«^(*2, *1) = S a p{q) + G^(je 2 , *i), (10) 

where q = X2 — X1, R = x 2 — Si, and x\ is the mirror im¬ 
age of X\ with respect to the plane 2 = 0 (see Fig.EKb)). 
The Oseen tensor given by 


FIG. 2. (a) A dumbbell in the simulated box with Lx Lx D. 

Periodic boundary condition has imposed at x-y plane, (b) 
A x-z projection of the dumbbell in the electro-osmotic flow, 
(c) An enlarged illustration of the bead in the dumbbell. It 
is composed of a large number of monomeric units of the 
polymers. 


(b) is the nonlinear electro-osmotic mobility. The mo¬ 
bility increases and is saturated with increasing electric 
filed. The threshold electric field Eq is typically 10 3 V/m 
that is experimentally accessible. 


Sa/3 (q) 


1 ( 

8tt?7o V 9 



( 11 ) 


where q is the magnitude of q, and the second term in 
eq. (ITOjl is 

G^(® 2 , *1) = -S a0 (R) + z 2 ( 1 - 25p z )V 2 R S a p(R) 

-22i(l -26p z )S ax ,p(R), (12) 


where 


Sa/3,7(^1) — ^g7So,^(q), (13) 


III. MODEL FOR SIMULATION 


and V g7 = d/dq^y. The Rotne-Prager approximation of 
the Blake tensor is given by^^r— 


In this section, our method of Brownian dynamic sim¬ 
ulation is described. As shown in Fig.[2](a), a dumbbell is 
simulated in a electrolyte solution with a no-slip bound¬ 
ary at 2 = 0. The dumbbell shows a dilute solution 
behavior. The solvent is described by a continuum fluid 
with the viscosity 770 and fills up the upper half of space 
(2 > 0). Electrolytes are also treated implicitly with the 
Debye length A = k - 1 . The dumbbell has two beads 
whose hydrodynamic radii are a, and actually the bead 
consists of many monomeric units of the polymer (see 
Fig. Etc)). The positions of the beads are represented 
by X\ and X 2 (see Fig. ED- Then we solve overdamped 
Langevin equations^ given by 

= uo(z n )S ax + J2 (G l^Fmp + k B TV mf ) G^) + 

m,/3 

for n = 1,2, a = x, y{&) 

where X na is the a component of a vector X n . uo(z) is 
the external plug flow, V na = d/dx na , G is the mobility 
tensor, F n = —V„17 is the force exerted on the nth bead, 

U is the interaction energy given as a function of x n . and 
k B T is the thermal energy. is the thermal noise which 
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1 + TrV? + -zr^l G°fl(® 2 , *1) + 0(a 4 ) 


for q > 2a, 


n RPB 

'-' r a/3 


(X 2 , X\) = < 


1 

6nrioa 



_9 q_(, 

32a ^ Q/3 


q a qp \ 

3q 2 ) 


+ ^1 + -g-V* + ^p(x 2l xi) + 0(a 4 ) 

for q < 2a. 

(14) 

The mobility tensor for the self part is given by£>2£ — 

G s „f (2) = lim G%0 B (x,xi) 

( M ||(*0 0 0 \ 

= 0 A*|| (z) 0 , (15) 

V 0 0 M-l(z) / 


where 


M|| (2) 
MJ l(z) 


1 

6irq 0 a 

1 

6nr/oa 


y a 1 /a\ 3 
~~ I62 + 8 \z) 
9a 1 /a\ 3 


1 - 


ya I /ay 
82 2 V 2 / 


+ ^(o 4 )) 

+ 0(a 4 ). 


(16) 

(17) 






























4 


Finally we obtain the mobility tensor as 

G™ = 6 nm G s $(z n ) + (1 - (18) 

The non-uniform mobility term in eq. © can be simpli¬ 
fied within using the Rotne-Prager approximation of the 
Blake tensor because the relation 

VmpG^ B {x n ,x m ) = 0, (19) 

fi=x,y,z 


is hold. Thus, the non-uniform mobility term is rewritten 
by 


55 ( 20 ) 

m,/3 

The interaction energy includes spring and bead-wall 
interaction given by 


U = U s (q) + ^ t/ w ( 2n ), 


( 21 ) 


n—1,2 

where U s is the spring energy as 


H 


U s (q ) = 


H 2, 

~~2 q ° ln 


-q , 


l - 


Hookian dumbbells, 


FENE dumbbells, 


( 22 ) 

where a FENE dumbbell stands for a finitely extensi¬ 
ble nonlinear elastic dumbbell, and a parameter b = 
Hq^/k-aT is defined for convenience. U w is the bead-wall 
interaction, 39 which is purely repulsive as 


U w {z) = 


2/a\ 10 /a\ 4 

5 W ~~ \z) + 


for z < a, 


(23) 


for z > a. 


Eq. © is numerically solved. Reflection boundary con¬ 
dition is set at z = D. When the center of the dumbbell 
goes across the boundary, the 2 -coordinate of each beads 
are transformed from z to 2D — z. For the lateral di¬ 
rections, the periodic boundary conditions are imposed. 
The size of the lateral directions is L x L. 


IV. RESULTS OF SIMULATION 

The concentration and velocity profiles are calculated 
by 



TABLE I. Simulation parameters. Nt is the number of total 
steps, Ni is the number of interval steps for observation, ]V m 
is the number of the sampling for each parameter set, and At 
is the time increment. 



Hookian FENE b = 600 

FENE b = 50 

N t 

5 x 10 10 

5 x 10 10 

25 x 10 10 

N 

5 x 10 3 

5 x 10 3 

25 x 10 3 

N m 

3 

3 

5 

At 

0.01 

0.0025 

0.0001 


where S(z ) is the delta function, Su(z) = u(z) — uo(z) is 
the velocity increment due to the polymeric stress, and 
(• • •) means a statistical average in steady states. The 
derivation of eq. (125|) is written in Appendix A. 

The imposed electro-osmotic flow uq(z) is given by 

U 0 (z) = MO E (1 - e ~ KS ), (26) 

where Mo is the electro-osmotic mobility in the pure sol¬ 
vent, and E is the applied electric fields Eq. m is 
rewritten in a dimensionless form with the length scale 
<5 0 = ^/k^T/H and time scale r = 6nr] 0 a/AH. The differ¬ 
ent types of dumbbells are simulated with the parameters 
noted in Table [I] It should be noted that the simulated 
systems are treated as dilute systems and the linearity 
with respect to the bulk polymer concentration is pre¬ 
served. After sample averaging, we obtain the concentra¬ 
tion at the upper boundary c(D), which slightly deviates 
from (A 2 !?) -1 because of the inhomogeneity near the sur¬ 
face. Hereafter, we define the normalized concentration 
as, 


C{z) 


c(-) 

c{DY 


(27) 


As well as the concentration profile, the velocity in¬ 
crement Su(z) has the linearity with respect to c(D). 
For convenience, we set a characteristic concentration 
Cb = O.l^o -3 , and the nonlinear electro-osmotic mobil¬ 
ity is defined by 


fi{E) = mo + 


Cb Su(D) 
c(D) E 


(28) 


The top boundary is placed at D = 100<5o, the lateral 
size is set to L = 1000<5o, and the Debye length is set 
to A = « -1 = 10<5o- We also set w = SfceT, and a 
hydrodynamic parameter h* as^i 


h* = —= 0.25. (29) 

V^oo 

Fig. [3] shows the steady state profiles of the Hookian- 
dumbbell concentration as functions of the distance from 
the wall. In the equilibrium state of E = 0, the profile 
shows the depletion layer whose width is of the same or¬ 
der as the gyration length <5 q - When the applied electric 
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FIG. 3. The concentration profiles of the Hookian dumbbells 
with varied applied electric fields. The inset shows the deple¬ 
tion length as a function of the applied field. The points are 
obtained by the Brownian dynamics simulation and the line 
is fitted by 8/5o = A(E / Eq) b , where A = 7.08 and B = 0.22. 


field is increased stronger, the depletion layer becomes 
larger than the equilibrium one and a peak is formed. 
The inset in Fig. [3] shows the depletion length as a func¬ 
tion of the applied field. The depletion length is defined 
by the position of the concentration peak. It shows a 
power-law behavior and its exponent is 0.22, which is 
much smaller than 2.0 in the case of a uniform shear 
flow^i The value of the concentration at the peak also 
increases as the electric field is enlarged. 

The results mentioned above are for the Hookian 
dumbbell which is infinitely extensible with the shear 
deformation. To consider more realistic polymers, the 
finitely extensible nonlinear elastic dumbbell is simu¬ 
lated. Fig. 0] (a) shows the concentration profiles at 
E = 1000I?o- Interestingly, the one-peak behaviors are 
also observed in the FENE dumbbells. In the case of 
the Hookian dumbbell, the concentration near the sur¬ 
face remains finite. On the other hand, in the case of the 
FENE dumbbells, the concentrations near the surface are 
negligibly small. Fig. 0 ] (b) plots the electro-osmotic mo¬ 
bilities with respect to the applied electric field. It is 
clearly shown the resultant electro-osmosis grows non¬ 
linear ly with respect to applied electric field. When the 
applied field gets stronger, the mobility increases and is 
saturated similarly to that in the toy model. The two 
types of the dumbbells have different rheological proper¬ 
ties from each other at the bulker— so that this non¬ 
linearity is not owing to the rheological properties of the 
dumbbells. On the other hand, the mobility is almost 
constant for E < lOEo, and this threshold of the linear¬ 
ity is larger than E 0 , that is predicted by the toy model. 
Likewise the saturation is observed when E ~ 10 4 Eo, 
which is larger than Ei. 

To clarify the difference of the profiles near the sur¬ 
face, ( q 2 ) and {q 2 /q 2 ) are plotted with respect to the 
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FIG. 4. (a) The concentration profiles of the polymers at 

E = 1000-Eo in the different types of dumbbells, (b) Plots of 
the nonlinear electro-osmotic mobilities with respect to E. 




(a) z / s o (b) z/s 0 


FIG. 5. (a) The profiles of (qt/q 2 ) at E — lOOOUo in the 

different types of the dumbbells, (b) The profiles of {(q/S o) 2 } 
at E = 1000i?o in the different types of the dumbbells. The 
curved lines are calculated by eqs. (IB 11) and (1B10I) . 


distance from the surface. Fig. 0 ] (a) shows the profiles 
of {q 2 /q 2 ). In the bulk, they approach to 1/3, which 
means the dumbbells are distributed isotropically. Near 
the surface, the polymers are inclined by the shear flow. 
Concerning the angles between the z-axis and the dumb¬ 
bell direction, that of the Hookian dumbbell is the largest 
among them. Fig. 0 ] (b) plots the profiles of (q 2 ). In the 
bulk, they approach to 3do which is the equilibrium value 
of them. Near the surface, they become larger since the 
polymers are elongated by the shear flow. In the case of 
FENE dumbbells, the saturations of the dumbbell length 
are observed. These behaviors are largely different from 
the minor difference in the concentration profiles. 


V. KINETIC THEORY 

In this section, a kinetic theory for a dumbbell is de¬ 
veloped based on Ma-Graham theory^. The probability 
function ^(aq, aq, t) obeys the continuity equation 

— = -Vi ■ (iq'F) - V 2 ■ (30) 

where x n is the flux velocity being given by^ 
x na = u 0 (z n )8 xa - J2 G ™ V ™p( U + ksTln^f). (31) 

m,(3 
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In the kinetic model, the beads are treated as point-like 
particles. Thus the mobility tensor is obtained by using 
G b instead of G RPB for both the self and distinct parts. 
The continuity equation can be rewritten with q and r 
as 

ST 

_ = -V r • (r*) - V 9 • ($*), (32) 

where r = (xi + a; 2 )/2 is the center of the mass of the 
dumbbell. We also define Vi and V 2 as 


Vi — — V r — V q, 
V 2 = — V r + V Q . 


(33) 

(34) 


Then, the probability function is also regarded as a func¬ 
tion of r and q. Here we neglect the interaction between 
the wall and beads. The flux velocities for r and q are 
obtained by 

fa = x[^o(£l) T ^0 {Z2^)]dxa "b ~zGapFp 


knT„ 


+ -^-G„^V^ln^-D^V rifl ln$, (35) 

qa — [^ 0 (^- 2 ) ‘Uo(zl')] 5 X a G a pFp 
krtT - 

+ ^-G a pV r ,f> In 4- - kBTG aP V qiP In T, (36) 

where F s = — Vi U s is the spring force, and D K is the 
Kirkwood diffusion tensor which characterizes the diffu- 
sivity of macromolecules, given by 


D 1 


fc B T(G n + G 12 + G 21 + G 22 ) 


(37) 


G and G are a variation of the mobility tensors defined 
by 

G = G 11 - G 12 + G 21 - G 22 , (38) 

G = G 11 — G 12 — G 21 + G 22 . (39) 

The concentration field c(r,f) can be obtained by in¬ 
tegrating the probability function with respect to the 
spring coordinate. It is given by 


c(r,t) = j 'Sf(r,q,t)dq. 


(40) 


We also define the probability function only for the spring 
coordinate as 




'&(r 1 q,t) 
c(r , t) 


(41) 


These new fields satisfy the continuity conditions, such 
that 


| = -V r • Mr),), 


(42) 


9T 

~dt 


= -V, ■(**), 


where (• • • ) q means the average with the spring coordi¬ 
nate, as 


(•■■),= 


/ ■ ■ ■'5>(r,q,t)dq 

c(r,t) 


■4'( r,q,t)dq. (44) 


For the limit of q <tC r, G and D K can be expanded with 
r. With keeping only the leading term, we obtain 


G = 


32nqoZ 2 


* ~q-z 0 ~XQx 

0 -g z ~xq y | + • • • > ( 45 ) 

V XQx XQy -2 qz 


d k = 


k B T 
12nqo a 


T 3 a _, , 

i + T s (q) 


where 


X = 


4 . ll+fi 

Az 2 


-5/2 


(46) 


(47) 


It should be noted that the approximation is more accu¬ 
rate than that in a previous studjdi since they considered 
only x ~ 1, which is not satisfied near the surface. With 
the approximation, ea. flMl) is averaged by T, and finally 
we obtain the concentration flux for 2 direction as 


c(fz)q — CU m ig(z) [ C (bK z )q] , 


(48) 


where 


nmigG) — k B TX7 q pG z p)q 


647T770 z 2 

x (x(q*F s x + q y Fy) - 2 q z F s z - 2 k B T( X - l)) q . 

(49) 

Eq. (gHD indicates two opposite fluxes of the polymers 
due to the external flow field. One is the migration flux 
from the wall toward the bulk and originates from the 
hydrodynamic interaction between the wall and the force 
dipolesdi The other is the diffusion flux from the bulk 
to the surface wall and is not found in the case of poly¬ 
mers in uniform shear flows dl It should be noted that 
the second flux includes not only the ordinary diffusion 
flux (D^ z ) q S7 r ,zC, but also the diffusion flux due to the 
q-inhomogeneity, cV r , z (D^ z ) q . When the external shear 
flow is uniform, the second flux vanishes, and the de¬ 
pletion length is proportional to the square of the shear 
rate since the migration velocity is proportional to the 
normal stress differenced! In the case of a plug flow, the 
diffusion flux suppresses the growth of the depletion layer 
and it may answer why the exponent of the depletion 
length is much smaller than 2.0 in the uniform shear flow. 
In steady states of the electro-osmosis, the total flux in 
eq. (1481) becomes zero, and thus, 


dc 

dz 


< Df z ) q 


'Urnie: 


d(D? z ) q 


(43) 


dz 


(50) 
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This equation shows the migration flux and the diffusion 
flux are balanced at the peak of the concentration pro¬ 
files. Finally the concentration profile is calculated by 


c = Cb exp 


Id (D*) q 



d(D? z ) g \ 
dz J 



The resultant flow profile can be calculated by 

5u{z) = f (z')dz', (52) 

Vo Jo 

where <x p is the polymeric part of the stress tensor as 


<r p = c(qF s ) q - ck B T I. (53) 

To obtain explicit expressions of c and Su , it it neces¬ 
sary to estimate u m i g , {D^ z ) q , and <r p . For this purpose, 
eq. 631) should be analyzed. However, eq. 61 is highly 
complicated. Even without the wall effects, it cannot be 
solved exactly, so that several approximation methods 
have been proposed* 44 ’ For simplicity, all the hydrody¬ 
namic interactions are ignored, and thus, the continuity 
equation is given 


(54) 

For the Hookian dumbbell eq. (l54l) can be solved for the 
second moment of q, and for the FENE dumbbell pre¬ 
averaged approximationdi’ 4 ^ is employed. The curved 
lines in Fig. [5] are calculated with these approximations, 
and they agree quantitatively well with the simulation 
results. In Appendix B, approximated expressions for 
these quantities of the Hookian and FENE dumbbells 
are written. 

Fig. [6] (a), (c) and (e) show the concentration pro¬ 
files for the applied field E = lOOOEo- The points are 
obtained by the Brownian dynamics simulation and the 
curved lines are obtained by the kinetic theory. The the¬ 
oretical calculations quantitatively cover well the simu¬ 
lations. Moreover, they reproduce the differences in the 
concentration near the surface between the Hookian and 
FENE dumbbells, since the migration velocities can be 
approximately proportional to (x)q (see Appendix iBl) . 
and it is much suppressed in the case of the Hookian 
dumbbells. Fig. [6] (b), (d), and (f) show the nonlin¬ 
ear electro-osmotic mobilities with respect to the applied 
field. The theoretical curved lines also have an acceptable 
tendency with the simulation results. However, they are 
not so consistent with the simulation results in weak ap¬ 
plied electric fields since the equilibrium depletion layer 
is not considered in the kinetic theory. 


d'F _ 

~ qa 


du 0 , 2 FI 

, dzOza r 
dz bTTTjoa 


T- 


2 k B T 
6nr]oa 


VI. SUMMARY AND REMARKS 

With Brownian dynamics simulations, nonlinear be¬ 
haviors of electro-osmosis of dilute polymer solutions are 



FIG. 6. (a) The concentration profiles of the Hookian dumb¬ 

bell as a function of the distance from the surface. The points 
show the simulation results and the curved line is calculated 
by the kinetic theory, (b) The nonlinear electro-osmotic mo¬ 
bilities of the Hookian dumbbell as a function of a applied 
electric field, (c) and (d) Those for the FENE dumbbell with 
b = 600. (e) and (f) Those for the FENE dumbbell with 

6 = 50. 


studied. The simulation results agree with a toy-model 
and analytical calculations of a kinetic theory. The main 
results are summarized below. 

(i) Under an external plug flow, the polymer migrates 
toward the bulk. The concentration profile of the 
polymer shows a depletion layer and a single peak. 
The thickness of the depletion layer depends on the 
electric field. At the peak, the migration flux is 
balanced to the diffusion flux. 

(ii) The growth of depletion layer leads to increment 
and saturation in the electro-osmotic mobility. 
Qualitatively this behavior does not depend on the 
rheological properties of the dumbbells. 

(iii) Analytical calculation of the concentration and 
the nonlinear mobility by the kinetic theory is in 
agreement with the Brownian dynamics simulation. 






















The threshold of the electric field for the nonlin¬ 
ear growth and saturation of the mobility is much 
larger than the prediction of the toy model, since 
the diffusive flux suppresses the migration toward 
the bulk due to the inhomogeneous shear flow. 


Here the microscopic expression of the stress tensor is 
given by 

= ~ 2 Fnm,a%nm,f3^ nm ( x )i (A3) 

n m^n 


We conclude this study with the following remarks. 

(1) Nonlinear electro-osmosis with A <C do has already 
been observed experimentally.—®' 27 They reported 
the mobility is increased with increasing the electric 
field. However, the nonlinear electro-osmosis with 
X^> So has not been reported experimentally, and 
therefore, experimental verification of our findings 
is highly desired. 

(2) It would be a future problem whether the hydro- 
dynamic interaction between the polymers and the 
surface plays an important role in electro-osmosis 
in polymer solutions even though A <C So or not. In 
this case the elongation of the polymers is strongly 
inhomogeneous under the plug flow with a short 
Debye length, and thus more realistic chain models 
should be considered. 

(3) Addition of charged polymers into solutions can 
change the direction of the linear electro-osmotic 
Sowj£i 22 When a sufficiently strong electric field is 
applied to this system, the direction of the flow 
might recover its original one. It needs to be inves¬ 
tigated theoretically and experimentally. 
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where F nm is the force exerted on the n-th bead from the 
m-th bead and S^ m (x) is the symmetrized delta function 
given by 


= / ds5(x — SX n — (1 — s)x m ). 

Jo 


(A4) 


The symmetrized delta function is integrated in the lat¬ 
eral directions as 


fi s nm ( z ) = / dxdy8 s nm (x) = / ds6(z - sz n - (1 - s)z m ) 


9(z - Zm) - 0(z - Z n ) 


(A5) 


Zn Zm 

where 9(z n — z) = 1 — 6{z — z n ). Then we obtain 

(z - z n )9(z - Zn) - (z - z m )9{z - Z m ) 


[ dz'Sg Tn (z') 
Jo 


Zm Zn 

min( 2 , z n ) - min( 2 , z m ) 


(A6) 


°n ^m 


where min (z,z„) = z9(z) — (z — z n )9(z — z n ). Finally, 
the velocity increment is expressed by 

1 f z - 

Su{z) = — J-J y ] film, 1 (~n — z m) / dz 5 nm (z ) 

= L2 Y, F nm>1 [min( 2 :, z n ) - min (z, z m )\ 


Vo L 2 


y-mm(z,z n )F nA . 


(A7) 


Appendix B: Approximated expressions for kinetic theory 
1. a Hookian dumbbell 


Appendix A: Derivation of the velocity equation for 
Brownian dynamics simulation 


Eq. (l54l) can be rewritten in a closed form for the sec¬ 
ond moment of the spring coordinates in a steady state 
with an imposed plug flow. The solution is given by^ 2 


In this appendix, the derivation of eq. (OS) is explained. 
The velocity field induced by the polymer is given by 

5u(z) = -—f (J% z (z)dz , (Al) 

Vo Jo 


(qq)g = 



where 


and the polymeric part of the stress tensor is obtained 
by averaging those of the microscopic expression in the 
lateral directions as 


J a/3 


= J 2 J dxdydlg(x). 


duo 

dz 


rnyoEe KZ . 


Therefore, we have 


(Bl) 


(B2) 


(B3) 


(A2) 


(qF s ) q = H(qq ) q , 
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and the polymeric stress tensor is 

<x p = c(qF s ) - ck B TI 

{ 2(t> 2 0 <j> \ 

= ck B T 0 0 0. (B4) 

U 0 o/ 

The Kirkwood diffusion constant can be estimated by 
k B T r ‘3- /i / 


<D* > g = 


12irr]oa 

k B T 

127r?7oa 


3a l ( q z 

i + t(A i + 7 

3a(q 2 +q 2 ) 


4 (q 2 )3/2 


(B5) 


where the second term is split into the second order mo¬ 
ments, and thus, we obtain 


(D?*) q = 


k B T 


1 + 


3a 2(</> 2 +2) 


q 127r?7oa [ 4<5 (2cp 2 + 3) 3 / 2 _ 

It is differentiated with z as 

-,K ' 


d K \ _ 3a 4ncp 2 ((p 2 + 3) 

dP zz)q ~ 12tt Vo a46 (2(/. 2 + 3) 5 / 2 ' 


(B6) 


(B7) 


The migration velocity can be estimated using the split¬ 
ting approximation of the averages as 

^mig(^) = 2 + <%) - 2x) 


64:7tt]oz* 
3 k B T 


y' q 

2 , J2 


(x)q(Qx+9 V - 2 )<5 


32tt770^2 


3k B T 


32nqoz 2 vx q 


(x)q<t > 2 


(B8) 


where 


(x)g = 


1 + 


1 + 


2 , 2 \ “ 5 / 2 \ 

\ 


4z 2 
+ 1 


2z 2 


(B9) 


where 


ib = 6\ —-— sinh < — arcsinh 
w 54 3 


b(j) f 3 + b 


108 V 54 


-3/2' 


The polymer stress tensor is 

( 2ip 2 0 ip 
cr p = ck B T 0 0 0 

^00 

The Kirkwood diffusion constant is 


(B12) 


(B13) 


(Dl) g = 


k B T 

q 127r?7oa 


j 3a Icp 2{ip 2 + 2) 


46 V ip (2^2 + 3) 3 / 2 


(B14) 


and its derivative is 


kBT 3a 
dz zz q 127rr7oa 46 


X K\ 


' dip 4ip(ip 2 + 3) 

(<pdip \ 

2(^,2 + 2 ) 1 

(2^ + 3)5/2 

\ip d(p ) 

(2lp 2 + 3)3/2 _ 


(B15) 


where 


~T1 = 2\/——7“ cos h \ -arcsinh 
dcp V 54 13 


b<p f 3 + b 


108 


tyV (b + 3^ 3 

108/ l 54 


108 V 54 
1/2 


-3/2' 


Finally the migration velocity is obtained as 

3 k B T 


UmiQ 


32'Krjoz- 


:{x)gfp , 


where 


(X)q 


ipip 2 + l\ 0/2 
<t> 2z 2 


(B16) 


(B17) 


(BIB) 


2. a FENE dumbbell 


The second moment of the spring coordinate for a 
FENE dumbbell can be obtained by pre-averaged clo¬ 
sures of p-FENE model4i>^2, It is given by 


and 


( Q<l)q = 


k B T ip 
H 6 


( 1 + 2ip 2 0 ip 
0 1 0 
ip 0 1 


^ 1 + 2ip 2 0 ip 


(qF) q = k B T 


\ 


0 

V> 


1 0 
0 1 


(B10) 


(Bll) 
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